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Most of the foundational work on quantum cosmology employs a simple minisuperspace model 
describing a Friedmann-Robertson- Walker universe containing a massive scalar field. We show that 
the classical limit of this model exhibits deterministic chaos and explore some of the consequences 
for the quantum theory. In particular, the breakdown of the WKB approximation calls into question 
many of the standard results in quantum cosmology. 
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It has been suggested by Hawking that the isotropy 
and homogeneity of the universe is a natural consequence 
of the Hartle-Hawking |^ "no-boundary" boundary con- 
ditions for the quantum state of the early universe. Test- 
ing this idea has proved difficult since a realistic model of 
the quantum state or wave function of the universe must 
describe the entire geometry and matter content of the 
universe. In order to make progress, a number of sim- 
ple models have been introduced. The most extensively 
studied model employs a homogeneous and isotropic 
Friedmann-Robertson- Walker (FRW) spacetime with a 
minimally coupled massive scalar field. Linde |^ showed 
that for certain initial conditions this model could lead 
to exponential expansion of the universe and a successful 
realisation of the inflationary paradigm. 

The FRW-scalar field model has been the work-horse of 
quantum cosmology. It is the model upon which most of 
the successes of quantum cosmology are based @-||], in- 
cluding Hawking's claim that the no-boundary proposal 
predicts inflation . In this letter we describe a problem 
that complicates the interpretation of these results. We 
show that the classical trajectories of the model exhibit a 
form of deterministic chaos known as chaotic scattering 
. This is a cause for concern since the quantum to clas- 
sical transition is qualitatively different for chaotic and 
non-chaotic systems. The problem is particularly pro- 
nounced in quantum cosmology since the theory is man- 
ifestly semi-classical and since the wavefunction of 
the universe is interpreted using the WKB approximation 
|,U|. As first pointed out by Einstein 0- the WKB 
approximation fails if the classical evolution is chaotic. 

Considering how much work has been done studying 
such a simple model, it may seem surprising that the 
chaotic behaviour went unnoticed. In fact, the finger- 
prints of chaos can be found in some of the early litera- 
ture pd|-p^, but it is only now, with an understanding of 
chaotic scattering, that we can recognise them as such. 

The model in question describes a closed FRW universe 
with Lorentzian metric 
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minimally coupled to a scalar field 4>. A fuller description 
of the model can be found in Ref. [^. The minisuper- 
space action is given hy S = ^ J {NH/a) dt where 



or' + m'^a^dP' is 



is the classical Haniiltonian and V - 
the minisuperspace potential. The classical equations of 
motion that follow from varying (p and a are: 
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Here an overdot denotes N~^d/dt. Variation of the lapse 
N leads to the constraint H = 0. 

The quantum description employs the "wavefunction 
of the universe" ^'(a, 0), which obeys the Wheeler- 
DeWitt equation = 0. Adopting a particular fac- 
tor ordering and changing coordinates to u = ae^'^ and 
V — ae"^, we find 
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where the minisuperspace potential now reads 

V{u, v) = uv ^-^uv {ln{v/u))'^ — 1^ . (6) 

In order to solve these equations we need boundary 
conditions. Hartle and Hawking have proposed that, in 
the Euclidean regime, the universe does not have any 
boundaries in space or time. The Euclidean version of 
the Lorentzian metric can be expressed in the form 



ds^ = dr^ 
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The no-boundary proposal demands that the geometry 
be compact and the matter fields regular. These con- 
ditions restrict the boundary conditions at zero volume 
(r = 0) to be 



0, 0: 
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By integrating the Euclidean action from t = to 
a nearby point (a, (f>) for large 0o a-nd applying WKB 
matching techniques. Hawking ^ obtained the bound- 
ary condition \l/ = 1 along the characteristics u = and 
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V = 0. We will subsequently question the validity of this 
procedure. The corresponding classical solutions take the 
form 



(m202)-i/2sin((m2</,2)V2^)^ ^ ^ 
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with (jjQ a constant. The Euclidean solution is matched 
onto the Lorentzian metric (|l]) near a totally geodesic 
spacelike surface S = 95"* = S'^ |Q where da/dr ~ 
0. The rotation to imaginary time is thus accomplished 
by taking r — -k / {2m(j)) + it, leading to the Lorentzian 
boundary conditions 
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Typical Lorentzian solutions see the universe inflate, 
then enter into a dust filled state where (j> oscillates 
rapidly, before recoUapsing to a final singularity as ^ — > 
±00 monotonically. In addition to these typical solutions 
there are an infinite collection of anomalous solutions [Q . 
The anomalous solutions expand, then recollapse, then 
bounce and expand again. The cycle of expansion and 
collapse may continue indefinitely, or may terminate af- 
ter a finite number of bounces. One of these bouncing 
solutions is shown in Fig. 1. Also marked are the lines 
V Q where the bounces occur, and the line uv = 1. 
It is difficult to track bouncing solutions for more than 
three of four bounces as they are exquisitely unstable. 
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FIG. 1. A bouncing classical trajectory in the {x,y) plane. 
Also shown are the lines V{u, u) = (dashed) and uv = 1 
(dotted). Here x = {v — u) /2 and y = (v + ii)/2. 

In a remarkable paper. Page |ll| conjectured that the 
anomalous solutions comprise a fractal set of perpetually 
bouncing universes. What Page described is now known 
as a strange saddle or strange repellor - the analog of 
a strange attractor for non-dissipative chaotic systems. 

The unstable periodic orbits partition the space of ini- 
tial conditions (here just the value of (/)o) according to 
their outcome. Unless a universe bounces perpetually it 



must eventually encounter one of the cosmological sin- 
gularities at a = 0, <?!) = ±00, corresponding to the lines 
u = 0, w = respectively. We can confirm Page's con- 
jecture by studying the boundaries between the two at- 
tractors m = and u = 0. We do this by color coding 
00 according to the outcome; grey for u — Q and black 
for V = 0. The attractor basin boundaries are shown in 
Fig. 2 for a representative range of (pQ. Each successive 
strip of "Universe DNA" resolves a small portion of the 
previous strip. We were able to track the fractal struc- 
ture over 12 decades in magnification before saturating 
machine precision. 
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FIG. 2. "Universe DNA": Zooming in on a boundary. 

In the close up views we start to see a third, white out- 
come occuring. These regions correspond to trajectories 
that inflate by at least 10 e-folds to produce a "big" uni- 
verse. A higher cut-off, such as the cosmologically inter- 
esting value of 60 e-folds, would produce much the same 
pictures, but with a smaller white region. The lower cut- 
off was chosen to keep the numerical integration short. 
It is interesting that trajectories near the bouncing so- 
lutions can lead to successful inflation (in the sense of 
solving the age and flatness problems) without requiring 
a large initial value for (j). For our choice of m = 1, a uni- 
verse that never bounces requires 0o > 5.88 to success- 
fully infiate, while a universe that bounces just once can 
successfully inflate \i 4>o ~ 1.19. Although the fractal set 
of perpetually bouncing solutions are countable and thus 
of zero measure, the collection of universes that bounce 
at least once is uncountably infinite and of non-zero mea- 
sure. For example, if we randomly chose an initial value 
of from the region shown in the upper strip of Fig. 2, 
a "big" universe is formed about once in every thousand 
attempts. 

The set of unstable periodic solutions can be quanti- 
fied by its spectrum of multifractal dimensions or by its 
topological entropy . Both of these methods provide 
a coordinate invariant measure of chaos in general rel- 
ativity pq. Here we calculate the topological entropy 
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as it can be found analytically. The topological entropy 
measures the growth in the number of periodic orbits as 
their period increases. We quantify the length of an orbit 
by the number of symbols needed to describe it. To do 
this we need to introduce a symbolic coding. The most 
efficient coding we could come up with records the sym- 
bol A for each upward crossing of the line uv — 1, and 
the symbol B for each crossing of the y-axis. For ex- 
ample, the bouncing trajectory shown in Fig. 1 has the 
coding ABABABBBABB. Applying this recipe to the 
first four primary orbits shown in Fig. 3 we obtain the 
codings 



1 = AB , 11 = ABB . 



Ill = ABBB , IV = ABBBB . 



(11) 



Here the overline denotes a sequence to be repeated. We 
can develop a recurrence relation for the number of pe- 
riodic orbits, N{k), with period k. Writing N{k) = 
P{k) + Q{k), where P{k) is the number of period k words 
ending in A and Q{k) the number ending in B, we find 

g(fc + l) = P(fc)+Q(fc), P(fc + l)=Q(fc), (12) 

and A^(2) = 2. The solution is then N{k) = (7*^+1 - 
'y^'^^^)/ \/5, where 7 = (1 + \/E)/2 is the golden mean. 
The uncountably infinite set of aperiodic orbits described 
by Page correspond to orbits with k = 00. 



I 




II 
























III ^ 




IV i 

















FIG. 3. The first four primary orbits in the {x, y) plane. 

Denoting the number of periodic orbits with lengths 
less than or equal to k by N{k), we find the strange 
repellor has a topological entropy of 



Ht = lim \ \iiM{k) = ln7 . 



(13) 



Since > 0, we can conclude that the dynamics is 
chaotic. It is interesting to note that the system has a 4- 
dimensional phase space - the minimum number required 
for chaos to occur in a Hamiltonian system. 



Having established that the classical evolution is 
chaotic, we can consider possible implications for the 
quantum theory ||l^,|l^. It has been suggested that the 
wave function of the universe will develop small scale 
structure if the classical dynamics is chaotic |l9|. This 
was demonstrated in Ref. by artificially choosing a 
potential with an infinite number of discontinuities. In 
the present more realistic model, we have to look else- 
where for the source of fine structure since the potential 
is completely smooth. Moreover, if we accept the bound- 
ary condition 5* = 1, there is no reason to expect small 
scale structure to develop since the evolution is governed 
by a linear wave equation with a smooth potential po|] . 

On the other hand, it is clear that the WKB approxi- 
mation breaks down for our model. This breakdown was 
independently observed by Shellard |l2) and Kiefer [|l3| . 
The WKB approximation seeks to interpret 5" in terms 
of classical trajectories according to the decomposition 



* = * I 



;C„exp(±i5„), (14) 



where each phase factor 5„ is taken to obey the 
Hamilton- Jacobi equation {W S)'^ + V = 0. The integral 
curves of VS* correspond to classical Lorentzian solutions: 

dS a . dS 

da 
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The WKB approximation is valid if the amplitude C 
varies much more slowly than the phase S. There is 
a conserved current J" = |C^|9"S'„ associated with the 
flux of each WKB wave packet. The fiux of a bundle of 
classical trajectories Fn = / J'^^a.pdx^ remains constant, 
regardless of the hypersurface on which it is sampled . 

In order for the flux to be constant, |C„p must in- 
crease if the bundle of trajectories focuses, and decrease 
if the bundle defocuses. When classical trajectories cross, 
a caustic occurs and |C„p 00. The basic WKB ap- 
proximation breaks down at a caustic, but an enhanced 
version continues to hold if the caustic is taken to be the 
source of a new complex conjugate pair of WKB solutions 
However, even this enhanced WKB approximation 
breaks down when the system is chaotic as regions near 
the unstable periodic orbits contain an infinite fractal set 
of caustics. The sum over the pairs spawned in these 
regions does not converge. In addition to the fractal set 
of caustics, the strange repellor also produces rapid de- 
focusing of nearby trajectories. This rapid divergence 
drives |C„p rapidly to zero, and the WKB approxima- 
tion again fails. Both of these effects can be seen in the 
numerical studies of Refs. p^ , p^ . 

At this juncture we appear to have a contradiction: so- 
lutions of the Wheeler-DeWitt equation with the bound- 
ary condition ^' = 1 appear completely benign, despite 
the classical evolution being chaotic and the WKB ap- 
proximation breaking down. 
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We believe the reason for the apparent contradiction 
can be traced to the boundary condition for In deriv- 
ing the condition = 1, Hawking only included a small 
subset of the non-singular bouncing universes. This was 
subsequently criticised by Page , who argued that the 
recollapsing, singular paths would modify Hawking's re- 
sult. Shellard jl^ attempted to calculate the contribu- 
tion from these paths using an enhanced WKB approxi- 
mation, but was thwarted by the caustics that developed 
near the unstable periodic orbits. 

Knowing that the root of the problem lies in the chaotic 
nature of the classical paths, we can suggest a new ap- 
proach. We conjecture that the contribution from the 
perpetually bouncing solutions, and the singular paths 
in their vicinity, can be found by performing a weighted 
sum over the multiple instantons that contribute to the 
bouncing solutions (see Fig. 4). The sum should be 
weighted by the instability exponent for each periodic 
orbit. The collection of Euclidean instantons that con- 
tribute to the strange repellor can be viewed as a fractal 
set of spherical "Russian Dolls" , with the radius of each 
4-sphere given by l/(r7i0o) where (j)Q takes all values on 
the fractal boundary. 
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FIG. 4. Some of the Euclidean instantons (solid circles) 
that contribute to the Lorentzian bouncing solution (dashed 
line) displayed in Fig. 1. 

Our proposal is motivated by Gutzwiller's ap- 
proach to quantum chaos. The Gutzwiller trace for- 
mula expresses the quantum-mechanical energy levels of 
a chaotic system in terms of a sum over the classical un- 
stable orbits. By adapting some of the techniques p^ ] 
developed to evaluate the Gutzwiller sum, we may be 
able to properly describe the fractal boundary conditions 
and subsequent chaos in the wave function. It would be 
interesting to study how these issues impact other ap- 
proaches to quantum cosmology, such as Vilcnkin's tun- 
neling proposal 1^. A breakdown in the WKB approxi- 
mation would complicate the notion of "out-going" wave- 



functions used in the tunneling approach. 

While our arguments are based on a particular model, 
in the context of a particular approach to quantum cos- 
mology, the issues we have raised will affect any the- 
ory of quantum gravity since all dynamical systems need 
boundary conditions, and generic dynamical systems are 
chaotic. 

We benefitted from discussions with Janna Levin and 
John Stewart. 
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